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Abstract. - The concept of an active zone in Laplacian transport is used to obtain an ana- 
lytical approximation for the reactive effectiveness of a pore with an arbitrary rough geometry. 
We show that this approximation is in very good agreement with direct numerical simulations 
performed over a wide range of diffusion-reaction conditions (i.e., with or without screening 
effects). In particular, we find that in most practical situations, the effect of roughness is to 
increase the intrinsic reaction rate by a geometrical factor, namely, the ratio between the real 
and the apparent surface area. We show that this simple geometrical information is sufficient to 
characterize the reactive effectiveness of a pore, in spite of the complex morphological features 
it might possess. 



The physics of diffusion-reaction phenomena in disordered media has attracted a lot of 
attention in recent years, due to its close connections with major technological subjects like 
catalysis and waste disposal. The development of modeling techniques for the description of 
heterogeneous catalytic processes represents a real challenge, mainly due to the limitations of 
classical pseudo- homogeneous representations These macroscopic models for diffusion- 

reaction processes can only implicitly account for the geometrical features of real pore spaces 

Due to screening effects, even if the active sites are uniformly distributed in space along the 
entire surface, the regions corresponding to fins or extended protrusions will display a higher 
activity as opposed to the deep parts of fjords which are more difficult to access. We recall that 
the standard modeling approach is to consider the catalyst particle as a homogeneous system 
where reagents and products can diffuse and react according to a given effective diffusion 
coefficient and an intrinsic reaction mechanism. Although in recent years some effort has 
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Fig. 1 - (a) Schematic representation of a two-dimensional rough pore. The reactive walls are com- 
posed of 10 successive different random Koch curves (RKC) of the third generation, (b) Schematic 
representation of the wall element. The reactive interface is irregular and the reactant concentration 
obeys the Laplace equation. A constant concentration Cs is maintained at the source line and the 
boundary condition at the interface is dC /dn = K S C / D. (c) Schematic representation of the coarse- 
graining or finite renormalization procedure used in the "Land Surveyor Approximation" (LSA). 
Successive ropes of length A determine successive chords of lengths L c i, L C 2,--. 

been dedicated to the investigation of diffusion and reaction in heterogeneous geometries 
(e.g., fractals) the role played by the pore shape or the local surface morphology on the 

overall diffusivity and reactivity of the catalyst remains to be understood. 

The notion of active zone || has been recently applied to elucidate several relevant features 
of the Laplacian transport to and across irregular interfaces (e.g., electrodes and membranes). 
In the present letter, we show that this concept can be used to derive an analytical approxi- 
mation for the diffusion-reaction activity of pores of arbitrary rough geometry. In a first step, 
we use a coarse-graining argument to obtain a general expression for the effective reaction rate 
of an irregular wall element. In the second step, we apply this result to calculate the reactive 
effectiveness of a deep rough pore analytically. Finally, this two-step approach is validated by 
direct numerical simulations of the full problem. 

As shown in Fig. [I], we consider the linear and steady state diffusion transport of reactant A 
in a single slit-shaped pore of nominal length Lt (x-direction) . The two reactive walls of this 
pore are rough surfaces of thickness b (z-direction) separated by a distance w (y-direction) . 
By definition, they arc generated through the translation along the z-direction of two lines 
made by the connection of several elementary irregular interfaces. In the example shown in 
Fig. [j](a) , the walls of the pore have a fractal geometry of a 3-generation random Koch curve 
(RKC) From a source at the pore entrance, molecules of species A diffuse into the pore 

through a neutral solvent to react at the walls according to a first-order kinetics, A — ► A* . 

The basic idea is to focus on a typical finite region of the pore enclosing a wall that 
corresponds to a single elementary unit of the interface. This subset of the system can be 
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approximately treated as a diffusion-reaction cell where reagent species A diffuse from a source 
line of length L towards an irregular reactive interface of perimeter L p and thickness b (see 
Fig. |f|(b)). In the bulk of this cell, the diffusion of species A obeys Fick's law, J(r) = —DVC, 
where J represents the mass flux vector field, G(r) is the local concentration at position r and 
D is the molecular diffusion coefficient. Under steady state conditions, the concentration field 
of the reactant species satisfies the stationary diffusion equation V 2 G = 0. For the boundary 
conditions, we assume a constant concentration of A equal to Cs at the source line. At the 
reactive interface, the depletion rate of A follows J n = —K S C, where K s is the intrinsic 
reaction rate constant. Due to mass conservation, the local reaction flux must be equal to 
the diffusive flux reaching the active interface from the bulk. As a consequence, the boundary 
condition at this point can be written as 



This introduces a finite length scale A = D/K s into the problem. Thus, in order to 
understand the diffusion-reaction behavior of irregular interfaces, one has to solve the Laplace 
equation with the boundary condition (Q) on the walls. This is an unsolved mathematical 
problem. 

In contrast, the problem of the Dirichlet boundary condition (G = 0) on irregular surfaces 
has been thoroughly studied, specially for the case of two-dimensional systems. An important 
theorem proposed by Makarov JllJ has been used to describe the properties of the current 
distribution on irregular electrodes (e.g, fractal electrodes) ||. This theorem states that the 
information dimension of the harmonic measure on a singly connected electrode in d = 2 is 
exactly equal to 1. In terms of activity, this means that, whatever the shape of the interface, 
the size L act of the region where most of the reaction takes place is of the order of the overall 
size (or diameter) L of the cell under a dilation transformation. When applied to a transport 
system like the diffusion- reaction cell shown in Fig. 0(b), Makarov's theorem has a simple 
but general consequence: the screening effect due to geometrical irregularities of the reactive 
interface can be characterized in terms of the ratio la] 



S = L p /L . (2) 

Because the active length of the interface is of the order of the system size, L act « L, then 
Lact ~ Lp/S and the factor S can be considered as the "screening factor" of the Dirichlet- 
Laplacian field. This result cannot be applied as such to a real diffusion-reaction cell because 
the boundary condition on the reacting interface is not C = 0, but rather dC/dn — G/A. To 
account for this more realistic boundary condition, we use the Land Surveyor Approximation 
(LSA) ||[l2|,|l3||. In this method, one switches from the real geometry obeying the real 
boundary condition to a coarse-grained geometry obeying the Dirichlet boundary condition, 
where the coarse-graining depends on A. The size of the active zone of the reactive interface 
can then be obtained from the screening factor of that new geometry. 

Consider a region i along the reactive interface with perimeter L pl = Sj+i — Sj, where s is 
the curvilinear coordinate along the interface. It has a reactive conductance Gi — [bK s }L p i. 
The access conductance to reach this surface is G acc w bD. Depending on the size of the 
region i, there exists two situations: Gi < G acc or G acc < Gi. If L pi is small, Gi < G acc and 
the flux is limited by Gi. If L pi is large enough, the flux is limited by the access conductance. 
In the latter situation we are, in a first approximation, back to the case of a Laplacian field 
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with Dirichlet boundary condition (i.e., C = 0). The new geometry is defined through the 
identity G; = G acc , which means that the curvilinear distance L p i along the reactive surface is 
equal to A. From its definition, a coarse-grained region can be considered as acting uniformly. 
If we performed the coarse-graining on a larger scale, the flux would no longer be uniform 
within a macrosite and, therefore, we would not be able to identify the size of the active zone. 
The size of a grain, L c , is the distance in real space (or the chord length) between si and 
S2- We then shift from the real geometry to the coarse-grained geometry which is made of 
successive chords L c i = L c (sq, s±), L c2 = L c (si, S2), ■ ■ ■ (see Fig. 0(c)). 

If we now observe that the perimeter of the coarse-grained interface is L p ^ cg — L c \ + L C 2 + 
... = N(L C ), where (L c ) is the average chord length, it is possible to express the screening factor 
of the coarse-grained geometry as S cg — N(L C )/L. In this way, the effective conductance of 
the interface can be written as G e s = bL p K s / S cg . By definition, since each chord corresponds 
to a perimeter A, we can write L p = NA and then express the conductance of the reactive 
interface as 



G c ff = bD . (3) 

\ L c) 

This result expresses how the conductance of a typical wall element should depend on the 
diffusivity of the reagent species in the bulk phase and the average chord length corresponding 
to a perimeter of length A. The geometry enters through the local relation between a perimeter 
of length A and its associated chord length. At this point, we can then define from @ an 
effective reaction rate constant for the interface as 



K cS = D/(L c ) . (4) 

This result is general. It applies to an arbitrary wall geometry provided that the wall itself 
does not contain deep pores [fl3|| . In the particular case where the wall irregularity can be 
described by a fractal geometry, the chord length (L c ) of a self-similar interface with dimension 
Df, size L and lower cut-off length scale i can be expressed approximately as 

r A if A < £, 

(L c ) = < i(A/£)^f if £ < A < L p , (5) 
[a/S ifA>L p . 

Note that in the case where A ^> L p , Eq. (^) becomes 



K cS = SK s . (6) 

This exact relation provides a very simple expression for the effect of the roughness on the 
effective reaction rate. 

In the second step, the LSA result (|]) is used to compute the overall activity of the 
rough pore by considering it as a smooth pore with an effective reactivity. Thus, writing 
this reactivity of the coarse-grained wall in terms of K e f[, it is possible to describe by mass 
conservation the phenomenon of diffusion and surface reaction over the entire pore using the 
following differential equation: 



J. S. Andrade Jr. , M. Filoche and B. Sapoval : ANALYTICAL APPROXIMATION... 5 



d 2 C d 2 C n 



with boundary conditions 



C(0,y) = C s , ^-(L T ,y)=0 , ^-( X> 0)=0 

and ^( X:W /2) + ^C(x,w/2)=0 . (8) 

The boundary condition ^(x,0) = corresponds to the assumption that the pore is 
symmetric with respect to y = 0. As in heterogeneous catalysis studies, the effect of reaction- 
diffusion process is measured by an effectiveness factor [^): 



where $ is the mass flux penetrating the system 



$ = 2b 



■j/2 



dC 



dy 



(10) 



and $ max = ZDbLxSCs / A is the mass flux in the absence of diffusion limitations. From 
the analytic solution of Eqs. (f?]) and (||) and using the definition (^), we obtain the following 
expression for the reactive effectiveness of a two-dimensional rough pore p4[ : 



2A ^ tan(a„it;/2) tanh(a„LT) 
71 = (L C )L T S ^ [ (a 2 + (L c y 2 ) W /2 + (L,)- 1 ] '"' 

where a n is the nth root of the eigenvalue equation, a n tan(a„w/2) = (L c ) 1 for n = 
1,2,3,... 

In order to check the validity of Eqs. (||) and ([ll]), we numerically calculate the solution 
of the two-dimensional Laplacian transport in a rough pore geometry whose active walls are 
composed by several elementary RKC interfaces (Df = log 4/ log 3). This task is performed 
here through discretization with unstructured meshes by means of the finite elements technique 
p5[. Due to computational limitations, we restrict our simulations to reactive pore walls made 
by RKC curves of 2 and 3 generations. In Figs. ||(a) and (b), we show the concentration fields 
in a typical reactive rough pore constructed with 3-generations RKC walls and calculated at 
A/i = 10 and 100, respectively. 

The results displayed in Fig. || show that the LSA predictions for the effectiveness are in 
excellent agreement with numerical simulations of the complete diffusion-reaction equations 
in both 2- and 3-generations RKC pores and over the whole range of relevant r\ variability, 
namely, for A > I jl(|. It is illustrative to compare the results obtained for a rough pore 
with the classical behavior of reactive, but flat pores. In the inset of Fig. ||, we show the 
logarithmic curves of r\ versus A/£ calculated for a rough pore (3-generations RKC) and for 
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Fig. 2 - Contour plot of the steady-state concentration of reactant penetrating a rough pore. The 
relevant sizes of the pore are w = 2L and Lt = 10i. In (a) A/£ = 10 while in (b) A/£ — 100. 
The contour plots have been computed from direct numerical simulation using the finite elements 
technique. The change in the grey scale corresponds to a factor of 2 in concentration. 

the corresponding smooth pores with lengths Lt and LtS. The large discrepancy which 
can be observed among the behaviors of rough and flat pores clearly indicates the important 
role played by the geometry of the reactive interfaces. Finally, one should note that, in the 
limit of negligible resistance to pore diffusion, the effectiveness approaches asymptotically its 
maximum, rj = 1. If the pore is very deep (Lt 3> u>), the one-dimensional approximation of 
Eq. (0) allows one to estimate the crossover for this saturation regime as A = L x = 2Lt 2 S/w. 
This prediction is also in good agreement with the numerical simulations (see the inset of 
Fig- 1). 
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Fig. 3 - Dependence of the effectiveness factor 77 on A/£. The symbols correspond to 77 values 
calculated from direct numerical simulation of reactive rough pores whose walls are made by RKC's 
of second (squares) and third (circles') generation. The thick solid lines correspond to the analytical 
approximations from Eqs. (^) and ([ll|). The arrows point to the crossovers A — L p and A = L x ■ For 
comparison, the inset shows two curves corresponding to the effectiveness behavior of smooth pores 
(dashed lines) with lengths Lt (top) and LtS (bottom). 
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In summary, the work presented here provides clear evidence that the effect of pore space 
morphology on the global efficiency of a diffusion-reaction system is not only relevant, but 
should govern the reactivity of irregular interfaces under diffusion limitations. In particular, 
it is found that, contrary to the classical description [jjj^], the effect of an arbitrary rough 
geometry is to modify the reaction rate and not the effective diffusion coefficient. Note that 
the analytical approximation derived here and validated through direct numerical simulations 
rely solely on the knowledge of simple geometrical characteristics of the interface to provide 
accurate predictions for the effectiveness of rough pores, over a wide range of diffusion- reaction 
conditions. For all practical situations where the reactant penetration in the pore is significant, 
one can use the simple model of a smooth pore with an effective reactivity K e g to describe 
the efficiency of a rough pore. Remarkably, K c g is the product of the intrinsic reactivity by 
the screening factor S, which has an elementary geometrical meaning. 

At last, since the Land Surveyor Method has also been successfully applied to the response 
of non- linear irregular electrodes [ fl7| , it could possibly be used to study higher order reactions 
in rough pores. 
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